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I. INTRODUCTION 


The study of the deformation of helical springs has most commonly been 
limited to those cases due to normal loading [Ref. 1, 4, 5, 6, 7, 8). In 
normal loading, a force is applied through the center of a spring, and the 
displacement of the spring is expressed as a function of the load and the 
parameters of the spring involving material properties and geometry. This 
thesis investigates the deformation of helical springs due to eccentric 
loading, with the aim of designing a helical spring that will deform to a 
given shape through the application of known forces. The spring design 
will include the material selection and the selection of helical spring 
geometry in terms of its length, the diameters of the spring and the 
spring coil, and the number of coils. 

Chapter II provides the background needed to study the deformation of 
the two helical spring models, such as the development of Castigliano’s 
theorem, and a brief introduction of Shape Memory Alloy’s (SMA’s). The 
descriptions of the two kinematic models that have been chosen for this 
study are also included in this chapter. Chapter III provides the 
analysis for spring deformation based on the first model, and Chapter IV 
provides the same for the second model. Chapter V contains computer 
programs that will be used for simulation based design of the proper 
helical spring. 

The motivation of this study is to design a helical spring that will 
be used as the primary structural element of a robotic manipulator for 
Minimally invasive surgical applications. In viewing the spring as a 
skeletal element, it is necessary to analyze the deformation of the 
skeleton under the action of external forces, so as to control its 
deformation or motion. The spring deformation will be produced through 


the use of Shape Memory Alloy (SMA) active elements which undergo a phase 


transformation upon heating above a certain temperature. When these SMA‘s 
undergo their phase transformation, they change their shape to a 
predetermined form. During this process they will deform the spring they 
are acting on. By properly selecting the arrangement of the SMA active 
elements about the spring, and controlling their phase transformation, a 


properly designed spring can be made to deform to a desired shape. 


II. PRELIMINARIES 


Certain preliminary topics need to be reviewed before an analysis of 
the deformation of helical springs due to eccentric loading can be carried 
out. This chapter includes the development of Castigliano’s theorem, an 
introductory examination of the properties of Shape Memory Alloy’s, and 


descriptions of the spring models we have chosen for analysis. 


AG CASTIGLIANO’S THEOREM 

Castigliano’s theorem is a very useful mathematical tool that can be 
used to study the deformation of elastic bodies under the application of 
generalized forces. The deformation of the elastic body is computed from 
the strain energy of the body. It is an energy based approach and therein 
lies its simplicity. 

Castigliano’s theorem, developed by Alberto Castigliano in 1879, is 
a method by which one can determine the deflection of an elastic body at 
the point of application of a force. If forces F, and F, are exerted on 
an elastic body at two different points, A and B, there are four 


associated deflections [Ref 2] 


Saat lak, (2.1) 
8 as=fasFs (2.2) 
Sna=faaF a (235) 
Sas=fask's (2.4) 


where the first subscript implies the point of interest, and the second 


implies the force of influence. The f's are constants, and are Known as 
influence coefficients. They represent the deflection of one point 
relative to the other. For example, f,, represents the deflection of point 
A relative to point B. These influence coefficients are properties of the 
elastic member. 


Maxwell’s law of reciprocity states [Ref. 2, p. 637] 


Fan=Fna (Ags) 


which implies that the deflection at point A due to a unit force applied 
at point B is equal to the deflection at point B due to the application of 
a unit force at A. 

Castigliano developed a method where the deflection due to multiple 
forces acting on an elastic body is obtained as the summation of 
deflections due to forces applied sequentially, one at a time. The final 
result is a set of equations like those found in (2.1) through (2.4). One 
must then superimpose these equations to obtain a series of equations for 
the displacements, or 6’s, at the different points where the forces act on 


the elastic body. For two forces F, and F, we obtain 


5 n=5 ant Sas =LaaF at laaks (2.6) 


6 5=5 pnt Onn= fay at lank s (2.7) 


Now one needs to determine the work done by each force at each point 


of the elastic body where it acts. The work done by F, at point A is 


Wra=—Fadaa (2.8) 


which, after substituting equation (2.1) becomes 


Wan=>Fa( fayFa) =4 faaFa (2-3) 


Likewise, the work done by F, at point B is 


Weems foes (2.10) 


If force F, continues acting on the body while F, 1s gradually 
applied, we see that there is additional work done on point A due to F,, 


namely 


Wan FO a5 i211) 


Similarly, if force F, continues acting on the body as F, is slowly 
applied, one can see that there is additional work done on point B due to 


F,, namely 


Waa=F'56 pa (2 aat2) 


After substituting equations (2.2) and (2.3) into equations (2.11) and 
(2.12), and by invoking Maxwell's law of reciprocity, equation (2.5), one 


obtains 


Was=lasF Fs (2716) 


Wea fapF aks (2.14) 


The total work done on the body, or its total strain energy, if F, is 


applied before F, is the summation of equations (2.9), (2.13), and (2.10): 


w O82 (ygPh +2 EyalaFgt Eas) (2.18) 


If force F, is applied first, and then F,, the order of addition becomes 


equations (2090), (2214), "anda 2) 


W=U=2 ( fag Fn +2 fag FF int lagFa) (2.16) 


Thus, the total work done on the body 1s irrespective of whether F, or F, 
is applied first. 
From Equation (2.16) one finds that the displacement at point A is 


equal to 


= aaFat fapF'n=5, (22155) 


and the displacement at point B is equal to 


OU =F ygFy+fpaFo=8, (2.18) 
B 


Castigliano’s theorem states that for any force F, acting on an elastic 
body, the deformation or deflection at the point of application of the 


LOrces wis 


- 9U 
8° Sp (2.19) 


in the direction of F,;, where U is the total strain energy of the elastic 


body under the application of forces. 


B. SHAPE MEMORY ALLOYS 

To cause the deformation of a helical spring one needs to apply an 
external force. Shape Memory Alloy (SMA) active elements were chosen to 
provide the necessary external forces. An SMA active element has low mass 


and a very high force to mass ratio; this attractive feature allows the 


miniaturization of the whole structure. A brief introduction of SMA’s is 
imperative to have a good understanding of how the spring will be 
deformed. 

A Shape Memory Alloy (SMA), is a metallic alloy that is given a 
certain predetermined shape at a high temperature. Once the alloy is 
cooled, it can be deformed, and will remain deformed until heated. Once 
heated above a certain temperature, the alloy “remembers” its undeformed 
shape and returns to it. 

There are many alloys that exhibit this shape memory effect. Among 
them are Ni-Ti, Ni-Ti-Cu, Cu-Al-Ni-Mn. The alloy is first shaped into its 
desired “undeformed" shape at a high temperature, when the microstructure 
is in its austenite phase. These “undeformed" shapes can vary greatly, 
but the primary shapes we are considering are those of a thin wire of a 
given length, or a rod with a given circular curve. Once formed, the 
alloy is quenched to allow the microstructure to come into its martensite 
phase. It 1s now ready to be deformed. 

The alloy can now be deformed by stretching it, bending it, or 
reshaping it by any one of a number of means. It will stay deformed from 
its original shape until it is once again heated up back into the 
austenite region, where it will return to its original shape. Te is 
beyond the scope of this thesis to present the microscopic analysis of 
this transformation. It is merely intended to explain what an SMA is and 
how it will be used as an actuator for spring deformation. 

The next two sections provide two spring models that describe the 
positioning of the SMA actuators relative to the spring. The SMA elements 
are in their deformed states initially, and they revert back to their 
undeformed states once heated above a certain temperature. During this 
process the spring is deformed. The undeformed shape or the memory shape 
can be given to the alloy by annealing for some time at a fixed 
temperature and then by rapid cooling back to room temperature. In the 


discussion to follow it will be assumed that the SMA has already been 


given its memory shape and attention will be focused on the design of the 


helical spring for achieving the design goal. 


GC. MODEL 1: HELICAL SPRING UNDER ECCENTRIC COMPRESSION 

This section considers a soring model consisting of a helical spring 
along with its actuators such that the spring can provide two rotational 
degrees of freedom besides a single degree of freedom for linear 
translation. In this model, shown in Figure 2.1, the helical spring is 
fitzred with two end caps. Three SMA wire actuators are attached to the 
end caps just outside of the spring. All of the SMA wires are to be of 
the same length, spaced 120° apart. During the process in which the SMA 
Wires are placed, the spring is given a small initial blas compression. 


This keeps the SMA wires taut and eliminates any slack in the wires. 





As current 1s applied to heat one of the SMA’s, the active element 
shrinks back to 1ts original "undeformed" length. During this process the 
ether two SMA wires remain in their deformed configuration, and the top 


plane of the spring bends over by virtue of eccentric compression. 


F,, F,;, and F, are forces exerted on the spring by the SMA wires, and 
F., F,, and F, are dummy forces. The deflection of the points of 
application of these forces can be readily obtained using Castigliano’s 
theorem. Position vectors, r,‘s, from an arbitrary point A to the points 
where the forces are applied are constructed. Angle 8 is a measurement 
taken from the point where F, is applied, around in a counter-clockwise 
manner. Using these position vectors, moments and torsions due to forces 
F, through F, are summed up at A. 

R is the radius of the spring, and L is the length of the spring. E 
in the modulus of elasticity of the spring material, and G is the shear 
modulus of that material. I is the area moment of inertia of the cross 
section of the spring coil, and J is the polar moment of inertia. The 
Humber Of spring coils is n. 

Using these values, the total stain energy of the spring can be 
calculated. Once this has been done, Castigliano’s theorem is invoked, 
and the displacement of the spring at any one of the six points of 
application of the forces can be found. Knowing the relative 
displacements of the different points on the spring coil, the angle of 
deflection can be computed. 

For the spring design problem, the angle of deflection is 
predetermined. When a helical spring is chosen, R, I, J, E, G, andn are 
known. From these quantities the force required to deflect the spring is 
calculated. If this force 1s one that the SMA wire can exert on the 
spring, the spring has been properly designed. If not, some of the 
parameters of the spring, geometry or material, must be changed and the 
forces recomputed. When the force required to deflect the spring matches 


the force the SMA can exert, the design problem is completed. 


D. MODEL 2: HELICAL SPRING UNDER BENDING 
In the second spring model, the assumption is that the spring bends 


under the action of SMA rods. Three SMA rods are attached to the spring 


internally along its length and placed 120° apart. Initially the SMAg@ea- 
ar2 in their deformed shape. When one of the SMA rods is heated it 


regains its undeformed shape and bends the spring in the process. 





Figure 2.2 Helical spring under bending. 


The SMA rod applies a force to each of the spring coils. Each of 
these forces are assumed to have two components, one normal to the coil 
directed towards the center of curvature of the spring denoted F, and the 
other in the tangential direction denoted f. From Figure 2.2, 98 is an 
angular measurement internal to the spring coil measured from the outer 
most point on the spring, moving in a counter-clockwise direction. R is 
the radius of the spring. The radius of curvature measured to the center 
of the spring is P, and p is the radius of curvature to the inside point 
of each of the coils as they are bent. Angle ® is the angle of curvature 
of halt of the spring. L is the length of the spring, nh is the numbereon 
coils, = is the modulus of elasticity of the spring material, G Teme 


shear modulus, I is the area moment of inertia, and J is the polar moment 


of inertia. For an arbitrary point A on the spring, r denotes the 
position vectors from A to the points of application of the forces, and @ 
is the angle subtended by each coil at the center of curvature of the 
spring. 

Consider now only the top half of the spring, since the spring is 
symmetric and the bottom half is identical to the top. When an SMA rod is 
bent through the application of current, each coil of the spring is acted 
upon by two forces, F and f. Using these forces, and the position vectors 
from point A to their point of application, a summation of all of the 
bending moments and torsions at A due to the action of these forces can be 
obtained. The total strain energy due to these moments and torsions is 
calculated, and Castigliano’s theorem is applied to compute the 
displacements at each point of application of the forces. 

To reiterate, for the helical spring design problem, the spring’s 
displacement is a given. Knowing the equations for the strain energy, and 
the material and geometry of the spring, one can work backwards to find 
the force required to bend the spring in this manner. One needs only to 
iterate using the spring material and geometry to find a_e force 


commensurate with the given displacement. 


ea 


III. DEFORMATION ANALYSIS OF THE FIRST HELICAL SPRING MODEL 


Recalling the description of the first spring model described in the 
preliminaries, one finds that F,, F;, and F; are forces exerted by the SMA 
wires on the spring. F,, F,;, and F; are dummy forces on the spring needed 
to find the displacements at their points of application. Point A is an 
arbitrary point on the spring, and @ is the angular measurement from F, 


arsuna in a countsr-clockwise —Girecc lon 


sete 
' 


COO | 


Figure 3.1 Helical spring under eccentric compression. 





There are s1x position vectors from the six points on the spring where 


the forces are applied to point A. 


F,=R(1-cos@) f-Rsin67Z (352)) 


N 


7,=-R(cos@-cos60°) 7+R(sin60° -sin6) 7 (33 


oF) 
ad 


T,=-R(cos60° +cos8) 7+R(sin6é0° -sinO) 7 ( 


F.=-R(1+cos6) i-Rsin67 (3.4) 


r,=-R(cos60° +cos6) 7-R(sin60°+sin®) 7 (3.5) 


7,=-R(cos0-cos60°) i-R(sinO+sin60°) 7 (356) 


Summation of moments about point A yields 


6 
Mant, 5+) (7,xF,) =0 (3.7) 
=1 
where, 
F,=-F,k P=) 2 (230) 


and where M and M, are the x and y components of the reaction moments at 
point A. 
After computing the cross products of (r, x F,), one can set the like 


vector components equal to each other and find 


M,=R[-(F,+F,) sin@+ (F,+F,) (sin60° -sin6) ae 
-(F,+F,) (sin60°+sin6) ] . 


M,=R(-F,(1-cos@) +F,(1+cos6) + (F,+F,) (cos@-cos60°) Geic) 
+(F,+F,) (cos6+cos60°) ] | 


To compute the strain energy, one needs to find the bending moment and 


torsion, m and €, on the spring coil. These can be obtained from the 


3 


reaction moments M, and M, through a coordinate transformation in the 


following manner: 


m\_ M,. (34a 
(2) =7(@)(¥ 
where 
r(0) = oon or ei 
-sin8é cos6 





Following the coordinate conversions, the moment and torsion equations 


become 


M=R((-F,+F,+ (-F,+F,+F,-F,) cos60°) sin6 alten 
+(F\+F,-F,-F,) sin60°cos6] €, 


C=R[ (F,+F,+F,+F,+F,+F,) +(-F,-F,+F,+F,) sin60°siné Bs 
+(-F,+F,+(-F,+F,+F,-F,) cos60°) cos6] €&, 


The shear force at any point of the spring coil can be obtained from 


the static equilibrium of forces: 


V= (F,+F,+F,+F,+F,+F,) K (3a 


The expression for strain energy due to shear is 
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L 
2 2 
u=f fv dx- f £0 rae (3.16) 
0 Q 


where U, 1s the energy due to shear, f is the form factor (10/9 for a solid 
circle), V is the shear force, L is the length of the spring (2nm), G is 
the shear modulus, A is the cross-sectional area of the spring wire, and 
x is the integration variable denoting the length of the spring coil. 
After making the substitution of J, the polar moment of inertia, for the 


circular area A by the expression 


and integrating with respect to @ from 0 to 2nm, the energy due to shear 
becomes 


Sl0nnr<-VeR aie) 


U 
¥ 18G7 


Snewcan now see that since r is much smaller that R, the r* term will 
dominate the numerator and make U, very small. Equation (3.17) will be 
compared later with the strain energy due to that of torsion to show this 
difference. 

The strain energy due to the bending moment 1s obtained through the 


equation 


6 2mm 
m? m? 

U_= ir ee C2443) 

«| ZEr J 2ETR® 


where U, 1s the strain energy due to the bending moment, m is the bending 
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moment, L is the length of the spring (2nt), and EI is the bending 
stiffness of the spring. After integration one obtains 


3 
y= TTR 


inp (Fr ty Fy-FyFy-2F, Fy-FyFo+F Fe +Fi+F,Fy-Fi 





-2F,F,-F,F,+F¢+FyF,-F,Fs-2FyFe+Fi (3.19) 


+F, F.-FyF,+Fs+FF,+Fe ) 


The strain energy due to torsion is obtained through the equation 


2nn 


L 
ST LCS ema. .o- 3-20) 





2JG 


where U, is the strain energy due to torsion, t is the torsion, L is the 
length of the spring, and JG is the torsional stiffness of the spring. 


Through integration the strain energy due to torsion becomes 


3 
Ue (3F24+5F, Fy +3F, Fy+2F,Fy+3 FF ot FF et3Fe 
+5F,F,+3 F,F,+2F,Fo+3 Fy Fet3 Fs +5 Fy Fy +3 FyF, (3.21) 


+2F,Fet3 Fe +5 Fi Fo+3FyFe +3 Fs +5 FeFe+3 Fe) 


A comparison of the strain energy due to shear and the strain energy 
due to the torsional moment U, demonstrates how much smaller the strain 
energy due to shear is. Assume that their are eight coils, a spring 
radius of 5 mm, a spring coil radius of 0.5 mm, that forces F, through F, 
are applied uniformly with a unit magnitude of 1 N, and that stainless 
steel is used as the spring’s material. The ratio of strain energy due to 


torsion, U,, to the strain energy due to shear, U,, is 


Vv 
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2 2 
= 4(4/; = poe sae a (seo) 


(F,+F,+~ +e)? 


Since this ratio 1s so large, one can see that strain energy due to shear 
1s very much smaller than the strain energy due to torsion. When the 
spring deforms by bending, the same can be shown for the ratio of the 
strain energy due to bending to the strain energy due to shear. This 
implies that the strain energy due to shear can be neglected. 

The total strain energy is simply the addition of equations (3.19) and 
ec), 


nnR?3 


U=U-+U,- SEI 





[Fi +F, F,-F,F,-2F,Fy-Fy FetFy Ft Fo +FoFy 


-F,F,-2F,F,-2F,F,+F3 +FyF,-FyF5-2F3Fe 
+Fi +F F.-F,F,+Fs +F5Fe+Fe | 


, ATR? (3.23) 


2 
Sog [SFI tS Fy Fo+3F Fy +2F, Fy+3 FF, 


+5 FF +3 Fp +5 F,F,+3F,Fy+2F,F,+3F Fe 
+3 F5+5F,F,(+3F,F,+2F,Fet3 FQ 


+5 FyF,+3FyFet3Fe+5 FFet3 Fe) 


By invoking Castigliano’s theorem the six displacements at the points 


of application of the forces are obtained as: 


icy 


Ne Uinta 


ET ar ae ae al 


PPE 


IG [6F, +5 F043 FF, +2F,43F +S FP, ) 


= nnR?* 


Spies Can oe Pcie chars | 


guriteles 


ae (OR tO Fa tS Fy +3 Fy +2F5+3 Fe] 


Dean ite 
> 2EI 





[-F,+F, +2 F,+F,-F,-2F,] 


pertain 


Soe (Ri tS Fa t6 Fy +5 PF. +3, +2 Fe) 


a fehiise 
of Er 





(-2F,-F,+F,+2F,+F,-F,] 


, mR? 


ae (AF t3 tS Fy +6 F,+5 Fo +3 Fe] 


eelidiag 


erie ey AI 


+ BEE 
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(3F,+2F,+3F,+5F,+6F,+5F.] 
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(3224) 


(229 


(3.226) 


(2200) 


(3.22S3) 





3 
6,= ae [F,-2F,-2F,-F,y+F,+2F,] 
(3 
3 
, OTR 


Say (SF t3F, +2 Fy+3 Fy +5 Fy +6 Fe) 





As stated before, F,, F,, and F, are dummy forces, which means 





























they do not really exist, they are used merely to compute 
displacements at those points. Setting F, = F, = Fe = 0 one obtains: 
8,= cc [2F,-F,- ~F,) + SEE [6 F,+3F,+3F,] (3 
§,= OUR (rp +F,- -2F,) + BEE (SF, +5 F,+2F,) (3. 
See : 
nk BRR (3 
8,= ae [=F 26. — Fe) + [pe t6e +35] 
_ nmR? ntR? 
Bee Wei a lst hs) + (2k, to8; O28: (3 
7 BERS [-F,-F,+2F,) + SAE. (37,137, +675] (3 
6 euti [F,-2F, oF) eee [SF,+2F,+5F,) (3 


Since all of the SMA actuators will be identical, the force they 


se) 


oo) 


that 


the 


ae 6), 


—32) 


ses) 


34) 


Lo 


will 


exert 1s expected to be the same. Consequently, i£ all of the Jemeee 
actuators are activated simultanecusly, the deformation at the different 
points on the spring will be the same. This implies that the spring will 
undergo compression. The interest here lies in the bending of the top 
plane of the spring. Therefore, all of the three SMA actuators will not 
be activated simultaneously. Now define Bh as the angle by which the top 
plane of the spring bends through when one or two of the SMA wires are 


activated simultaneously. 





Figure 3.2 Deflection angle of the top plane ot the spring, 6. 


The two possible cases to cause the bending of the top plane of the 
spring are the application of one force, and the application of two 
rorces. In the first case one force, F,, is applied. This single force 
will cause displacements to occur at all six points on the spring. Of 
interest are the deflections at points 1 and 4 to determine the bending 
angle of the top plane of the spring. The deflections at points 27a 
and 6 will be such that the bending will not take place at in any other 
d2reetre tin By knowing the displacements at points 1 and 4, one can 


calculate B. With F, applied one finds 


Figure 3.3 Angle of detlection, 
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(3 om) 





8-8 
kang 2s (3.233. 
which means 
Z a al 
=tan UnnrR?F,/ — +— (22) 
‘ ee (er sa 
When any one force, say F,, is applied on the spring, the angle B is 


given by 


to 
' 


oe 1 al 
B=tan [nmr r, Sawa (3.40) 


The other important case is when two forces are applied 


Simultaneously, say F, and F,. In this case the important displacements 


are at points 2 and 5, where 


= aL 5 
& =ntR7(F +F) 71,3. (3.42) 
: a i 2bie 206 


This implies that for the application of the two forces F, and F; 


B=tan”| > nmR? het) all (3.43) 


Since each of the SMA wire actuators will be identical, the force that 
each will exert on the spring will be equal in magnitude. This means, 


then, 1f two SMA wires are actuated the resultant angle of deflection will 


be 


Spal 2p{_2 ,_1 3.44 
B=tan [nn ae 5a) ( ) 


From these results one finds that if one SMA wire is actuated, or two 
are, the angle of deflection is the same. The difference is the direction 
in which the helical spring bends over. It is now clear that by 
activating one or two actuators it will be possible to bend the top coil 
of the spring in three separate directions by a positive or a negative 


angle B. The case of actuating all three wires has not been investigated 


22 


Since this would produce the same result as that of applying a single 
aXlal force through the center of the spring, or applying a normal load. 

Now that the effects of helical spring geometry and material 
properties are known in computing the deflection angle fB, use of these 
relationships will be instrumental in designing the proper helical spring. 
SMA wires deform by a constant amount of almost 5% of their original 
length in the presence or absence of external forces. However, for this 
property of the SMA to be exhibited repeatedly, the opposing stresses in 
the SMA wires should not exceed a certain value. The design of the 
helical spring will be computed from a known value of $B that will give the 
force required to produce that deflection. If these forces produce 
stresses that are below the stress limit of the SMA wires, then the design 
1s feasible. However, one should try to increase the stresses in the wire 
as high as possible, within limits, such that the spring to be designed is 
not too soft. A spring that is too soft will not be strong enough to act 


as a manipulator for surgical applications. 
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IV. DEFORMATION ANALYSIS OF THE SECOND HELICAL SPRING MODEL 


The second helical spring model considered is one where the spring is 
bent over directly by an SMA rod. Three rods will actually be placed 
inside the spring 120° apart to give the spring a full range of motion. 
Here, consideration 1s given to only one SMA rod for the development of 


che model. 





Figure 4.1 delical spring under bending. 


The only major assumption made is that the length of the spring is the 
same in the bent over configuration as in the underormed configuration. 
With this in mind, certain quantities are defined. Let R be the radius of 
che spring and L be the length of the undeformed spring. in the deformed 
configuration, it is desired that the spring bends into the shape of a 


Circular arc. Let P be the radius of curvatures of the central ax: s oa 


spring, and p the radius of curvature of the inner part of the spring. 
Let ® be the angle between the top coil of the spring and the central 
coil, 0 be the angular change of spring coils, @ be the angular 
measurement along the length of the coils of the spring, @ be the angle 
subtended by two adjacent spring coils, and n be the total number of 
coils. From these definitions, one sees from Figure 4.1 that the 


following identities hold: 


o-—, (can) 

a- 22 (4.2) 
Nn 

o-(> |e (4.3) 

p=P-R (4.4) 


In this model there are three different coordinate frames that are 
used. The traditional f, f, and R frame is fixed on the top coil of the 
spring. The f’', #’, and R’' frame is relative to each succeeding coil. 
The €, and €, frame is a normal and tangential transformation of the f, f, 
and R frame. 

Suppose that as the SMA bends over, it exerts two forces on each 
spring coil it touches, F in the normal direction, and f in the tangential 
direction. The line of action of the normal forces passes through the 
center of curvature. For the ease of computation, it is helpful to 
translate the point of application of the normal force along its line of 
action to the center of curvature of the spring. Point A, at the top 


center of the spring, is the point where position vectors from all other 


Pa) 


points on the spring will be constructed to. In considering only the 
first coil, there are two position vectors from™each point on thaegecwe 
one corresponding to the normal force F, and one to the tangential force 
f. These position vectors have been defined this way so that they will 


apply to each coil for these forces. They are 


7,=(P-Rcos®) i “-Rsin6j’ (4.5) 


7, ‘=(P-Rcos®) f’-Rsin®j*-pi (4.6) 


The primes on the 1 and j vectors indicate that they are in the spring 
coil’‘s frame of reference, where the unprimed vectors are in the fixed 
frame of the top coil. These vectors will give the position from any 


point on the first coil relative to point A. After recognizing that 


= 


F\= Fi 


F,cos@i “-F,sinok’ (4p 


f,= £,k = £,sinoi’+f,cosok’ (4.8) 


the total moment about point A caused by these forces can be calculated by 


M, = (2,xF,) +(7, “xf,) (4.9) 


For the second coil the position vectors remain the same, but the 


forces acting on this coil are F,, f,, F,, and f, where 


F,=F,i=F,cos (a+) i ’-F,sin (a+) k’ (4.10) 
F,=f,k=f£,sin(a+o) I°+f,cos (a+) k’ (4.11) 


Zo 


F,=F,cosoi “-F,sinok’ (4.12) 


F,=f,sinoi “+f,cosok’ (4.13) 


Now the total moment about point A is 


eS (ahh lon) ies, se eal Gelato ay), (4.14) 


By proceding in a like manner, one arrives at the k-th coil, where 
Mpeemaowleils acted upen by the forces: FP), Fo; slo; eFp; and fyofok 5.7 ie 


These forces take the form 


F,=F,cos[(k-1) a+] i°-F,sin[(k-1) a+] k’ (4.15) 
F,=F,cos [(k-2) a+] i°-F,sin[(k-2) a+] k’ (4.16) 
PF =F,cosoi “-F,sinok’ (4) 
fF,=f,sin[(k-1) a+] 1°+f,cos [(k-1) a+] k’ (4.18) 
£,=f,sin[(k-2) a+) f°+£,cos [(k-2) a+] k’ (4.19) 
£,=f,sinol “+f,cosok’ (4.20) 


Now the total moment about point A becomes 
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k k 
Ok, zor) ye eG) (4.21) 


or, after computing the cross products 


k 
De, {Rsin@ (F,sin[(k-1) a+] -f,cos[(k-1)a+o]) i’ 


+[(P-Rcos®) (F,sin[ (k-1) a+] (4.22) 
-f£,cos [ (kK-i) a+] ) +£,p) 7° 
+Rsin@(F,cos [ (k-i) a+] +f,sin[(k-i) a+] ) k*} 


The easiest way to work with this moment is to break it into parts as 


follows 


k 
2 {F, (-Rsin@sin [ (k-i) a+] ) 
=1 


(4.23) 
+f; (Rsin6cos [(k-i) a+] )} 
k 
aa, {F, (-(P-Rcos®) sin [(k-i) a+@] ) em 
+f, ( (P-Rsin@) cos [ (k-i) a+] -p)} 
> UF 
M,-= -Rsin@ k-i) a+ 
: > ,(-Rsin@cos [ (k-i) a+] ) “aoe 


+f, (-Rsin@sin [ (k-i) a+] )} 


Zo 


A coordinate transformation is now in order to transform the moments 
from the coil fixed frame to a normal-tangential frame. The 


transformation matrix T(@), from equation (3.16), will be used such that 


(2) (* | (4.26) 


cos8 sin#é 
-sin6é cos6 








The moment equations (4.23) through (4.25) become 


k 
o> {F, (-Psin@sin [ (k-i) «&+@] ) (4.27) 


+f, (Psin6cos [ (k-i) a+] -psin6)} 


k 
Me= 9 {F, ( (R-Pcos6) sin [(k-i) a+] ) (4.28) 


+f, (-(R-Pcos®) cos [ (k-i) a+] -psin6)} 


k 
ye Re ese) C220) 


+f,(-RsinOsin [ (k-i) a+] )) 


With each of these moments in the proper frame of reference, the 
strain energy due to each of them must be calculated. With the subscript 
k on each of the energies to remind one that this 1s the strain energy of 


the k-th coil, the strain energy due to each of these moments is given by 
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U, -=| ==> Rd (4.30) 
) 


U,. = [ —t.RdB (4.31) 





2n 
ay LB (4 22o5 


Since each of the moment terms are expressed as a sum of k quantities, 
the integrals involving the square of the moments becomes quite 
cumbersome. To simplify the computation, one needs to use matrix algebra. 


First define the incremental force vector as follows 
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AF, 
AF, 


AF, 

AF=| - e R2* (4233) 
Af, 
ie. 


Af, 


Then define 


-Psin@[ (k-1) a+] 
-Psin@sin [ (k- 2) a+] 


Ae « Rt (4.34) 


-PsinOsind 


Psin@cos [ (k- 2a -psiné 
Psin®cos [ (k-2) a+] -psin® 


A,= e R* (4. 35) 


Besnbcas sensane 


It can now be shown that 
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A, 
M_=F7| . (42316) 


Likewise for the torsional moment, define 


fe Re So aed 
(R-Pcos®8) sin[ (k-2) a+] 
B= : (43,295 


= 


(R-Pcos®8) sind 


- (R-Pcos®) cos [ (k-1) 4 eee 
- (R-Pcos®8) cos [ (k-2) a+] -pcos6 
B= (4.38) 


- (R-Pcos®8) cos$-pcos8 


Then, 


B, 
M,.=F 7 (4.32) 


For the moment in the z’ direction, define 
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ere oe eee 
-Rsin@Ocos [ (k-2) a+@] 
C,= (4.40) 


-RsinOcoso 


Be get aie 3 een 
-RsinOsin[ (k-2) a+] 


C,= (44 1) 
-RsinOsind 
Then, 
Cy 
M,-=F?| (4.42) 
Cr 


With the moments defined by Equations (4.36), (4.39), and (4.42), one 


can obtain their squares as follows: 


A, 
weeart{*}aiay ‘AF (4.43) 
7 


B, 
Me=AFT -. 





}8,)aF (4.44) 
2 
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Cy 
M-=AF DA ees 
C 


2 


|e ey AF 


Now define 


A,A, } A,A,| |A,. ¢ Ap 
°(A,:A,) = a P|) |e € R2*kx2k 
Aya, AyAd| An Abs 


A, 


eee 


A= 








one arrives at 


MZ=AFT-A-AF 


where the elements of A,, are of the form 


A,, (i,j) =P? sin?® sin[(k-i) e+) sin[ (k-j) a+) 


the elements of A,, are of the form 


A,,(41,j) =-P* sin’6 sin[ (k-1) a+] cos [ (k-j) a+] 
+Ppsin’6 sin[ (k-i) a+] 


Ao) At 


and the elements of A,, are of the form 
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(4.45) 


(4.46) 


(4.47) 


(4.48) 


(4.49) 


(4.550) 


A, (1,7) =P’ sin*6cos [ (k-1) a+] cos [ (k-j) a+] 
-Ppsin’@ (cos ([(k-1i) a+] +cos[(k-7) a+] ) (4.51) 


+p*sin’6 


Similarly, by defining 


B, B,B, : B,B,} |By, +} By 
wo 1e(Bit Ba) = ee PRP we eae (4.52) 
72 B,B, + BzBz) [B21 * Baz 


B= 





one can see that 


Mf=AFT BAF (4.53) 


where elements of B,, take the form 


B,, (i,j) =(R-Pcos®@) *sin [ (k-1) a+] sin[ (k-j) a+] (4.54) 


the elements of B,, take the form 


B,, (1,j) =- (R-Pcos®@) *sin [ (k-i) a+] cos [(k-j) a+] 
- (R-Pcos®8) pcos@Osin[ (k-i) a+] 


(42:55) 


B,, (i,j) =By, (4,5)? (4.56) 


and the elements of B,, take the form 
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B,, (4,7) = (R-Pcos®) cos [(k-1) a+] cos [(k-7) a+] 
+(R-Pcos®) pcos6 (cos [ (k-i) a+] +cos [ (k-j) a+] ) 


(4057) 
+p2cos*6 
By defining 
‘a CC, §O,C, Cir? Che 
of: (Cae = tee ae. —) rr (4.58) 
Ce C20, 'O,0, Coz § Coz 
one can show that 
M.-=AFT-CAF (qe 
where elements of C,, take the form 
C,, (i,j) =R*sin*Ocos [ (k-i) a+] cos [ (k-7) a+] (4.60) 
the elements of C,, take the form 
C,, (i,j) =R*sin’O6sin [ (k-i) a+] cos [ (k-7) a+] (4.61) 
Cp, (4,65) = Cy, (1,5)* (4.62) 
and the elements of C,, take the form 
C,, (i,j) =R*sin*Osin [ (k-i) a+] sin[ (k-j) a+] (4-635 


Now that the squared moments are in compact form, the integrals can 
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be evaluated for the computation of the strain energies. The strain 


energy in the k-th coil due to the bending moment M, is 


2n 


(ex 


M2 


: (4.64) 
na [fh ae 


Rdd= Ryo h 








where one must evaluate the integral of A term by term. This means that 


te sone aefines 


2K 
Aas) 6 {(gez}an 4D) de 





_ (4.65) 
- oh, [PPsin'esin  (K- i) a+] sin[(k-j) a+] ce 
_| PYR ae ae x we 
bey Peay sinacos [a (2k-i-j+1)]+ 5 COS [a (i-j)] 
2K 
Ay, (4,5) af (555 }Ar (4,9) de 
0 
Bigs 2n-r- : 2 Wests 
(ser) ocar-qnay Sina sin la (2k-3 j+i)] ates 
-F sina (i-J)] 


-— 328 Fe —ein( S)sin (a (k- -it+=]} 


a (a*-1677) 


oe 


A; Aig) 6A (AGa) 


and 


Pde 3 
Foe “(einai 6 (gex}ee (4-9) de 


2[ OR \;2 eens eee 
eal aw (a*-417) aa ae (4.68) 
+E cos [a (i-Jj)] +np? 


64n°P wn flee a Te ee 
ee eazy Stn{ 5 )e0s [-S (2k-4-F41) 1cos (5 (4-9) ) 


then, the energy in the k-th coil due to the moment M, is 


Ay, “3 A, ° 
U, APTA CAF=AFT™ ~~~ LAF (4.69) 
An,” Aa?” 


The energy in the k-th coil due to the torsional moment M, takes the 


form 
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2% 


jePo 


2n 2 


M 
U.. ko Bee Ar 
fe] 


oes (4.70) 








where a term by term integration, as before, yields Equation (4.75), where 


the elements of B,,’, Bi2’, Ba,’', and B.’ are defined as follows: 


B,, (i,j) “(58 5)'n(R2+E Joos (a (i-3)] 





CAL) 
TR? 2naRP, nP* (a*-227) . 
{22° - 2208? , xP (?-20) )ping COs [a (2k-1-J +1) }} 
PEA aR 2, FP). ae 
Beoti7 3) (2 )ta( ic |sinta(s J)] 
{22- 2T7aRP , TP? (a?-2707) sina sin [a@a(2k-1-j+1)] (4272) 
a a?-n2 2 = a (a*-4 277) 
4n(a?-8n*)Pp_ 4naRp \_:_/a\_. Facies 
[sete fe ae sin{ > )sin [a (k i+2)}} 
Bo (ta) Sos (4273) 


and 


22 


Bayo -(58,|'2{22+F Joos [a (i-j) ] +n? 


o( BE - 208P « 2P° (928°) Join cos [a (2k-i-j+1) ] 


a a*-n2 (a*-477) (4.74) 


{s2ie 877) Pp _ snake |sin(S)cos (a (x- -i+2)]} 


a (a*-1672) a*-4n? 


{sz tal-en)Pp_ 4nake \sin{ #)cos [a (k-j+4)]] 


a (a?-1677) a?-4n2 


Thus the strain energy of the k-th coil due to the torsional moment M, is 


By,” ? By” 
U, ,= AFB -AF=AF?: ies.) -saullAue (4.753 
Ba,” ¢ Ba” 


The strain energy of the k-th coil due to the moment M,. is 


—2" RdB=AFT: f( Egle} ar ‘AF (4.76) 


taf HE 








where the elements of the submatricies C,,’, C,.’, Cs,’, and C,,‘ are defined 


in the following manner: 
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sinacos [a (2kK-1-7+1) ) 


R? ». Zar 
a (a*-477) CAgn a) 


ea) 
see) e 


+ Ecos [a (i-j) }} 


ca (4.9) BE -— 2 sine sin ta (2k-4-9+2)1 
- Ze a (a*-477) (473) 
+Zsin[a(i-j)]} 
Grr =e ta) (4.79) 
and 
ae = 2%? ee 
Ge 1,3) = — |\— mS Sinacos (a (2k-1-jtl1) ] 
99 (4.7) a rea ers 3-1) ines 


+5cos [a (i-j)])} 


Thus the energy of the k-th coil due to the moment M,, is 
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C1” : Cy2° 
U,- ppAPT C AF=AFT eee sal We (4.81) 


The total energy of the k-th coil of the helical spring is the 


summation of the energies due to the three moments, or 


Gly 


n 


tO, .tU,- AFT (A +B°+C’) AF (4.82) 


If one defines 


DaA’+B°+C° (4.83) 


the incremental energy of the helical spring up through the k-th coil now 


becomes 


k k 
Au-y 4) (AFTD-AF) , (4.84) 


Castigliano’s theorem is now invoked to determine the incremental 
displacements of each of the coils in the normal and tangential 


adirections. 
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b( BB) =DysAF, +2D,,AF, +2D,,AF*. ete Py 





AU \_ 
8 Je2D,,AF, +D,.BF, +20, )AF +. Be OP OO 





7 BE ped iaA A +20 AF, +20 QAP, + .. +D AF, 


k 


which simplifies to 


8, D,, 2D,;z 2D,;,; °° * 2D, AF, 


oe _|2P2: Dz2_ 2D2; °° * 2Dy AF, 


: he 
21) Os RN) 9 aes a mE 2) : 


Now define 


D,, 2D,2 2D;; °° * 2D, 
Alea Doz 2D2,; °° * 2D, 


Di ie Dees eat 


then the final important relationship becomes 
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(4.85) 


(4.86) 


(4.87) 


(4.88) 


(4.89) 


AF=A7*8 (4.90) 


With this relationship, knowing the physical geometry and material 
properties of the helical spring, and specifing a desired displacement, 
one can determine the force required by the SMA active element to bend the 
spring. The matrix A contains a conglomeration of information based on 
all of these known quantities. Once these given quantities have been 
specified, it is a simple matter to determine the required force matrix F. 
If the force required cannot be achieved by the SMA, iteration 1s required 
for the spring’s geometry or material properties. A program is enclosed 


that will aid in the design of a helical spring under these conditions. 
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V. SUMMARY AND RECOMMENDATIONS 


This thesis investigated the deformation of helical springs under 
eccentric loading using two different models. The goal was to develop an 
algorithm for computing the deflection of a helical spring under a known 
force. The desired final shape of the spring was that of a spring bent 
over in a circular fashion. The two models used achieved this by two 
different methods, both utilizing SMA active elements in two different 
configurations. 

The development of this study was to obtain relationships between 
force and displacement so that one could predict or control the 
deformation of the spring. A spring under normal loading is governed by 


an equation such as 


F=k6 (Sly) 


where force F is proportional to displacement 6. This proportionality 
constant, k, 1s based on the spring geometry and material properties. In 
the development of this study, similar relationships were obtained, but 
with a more complex form of proportionality factor. This factor ls 
dependent on the spring configuration, and 1s not a constant. The 
proportionality factor is also a function of the geometry and material 
properties. 

This study has been a first step in the design of a robotic 
manipulator to be composed of a series of helical springs that utilize SMA 
actuators to control the manipulator. In the course of future research, 
SMA actuator design is needed, and in assembling the helical spring with 
the actuators. Once the actuators are designed, the results acquire by 


this study can be employed to match a spring’s geometry and material to 
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the force that can be generated by the SMA to bend the spring in a desired 
manner. 

When the design process is completed, a helical spring and actuator 
system will be realized, and can be implemented in a minimally invasive 
Surgical manipulator. The conclusions of this study can also be used in 
any other situation where an eccentric load is placed on a helical spring 


causing it to bend over. 
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PROGRAMS AND PLOTS 
The first algorithm delevoped was for the first model, where the 
spring was deformed by applying an eccentric load axially on the spring. 
It is a FORTRAN program where material properties, spring geometry, and 
desired deflection angle are required as inputs, and the force required to 


deform the spring is the output. 


PROGRAM THESIS 


* This program is for the spring design used in the manipulator skelton 
* for minimally invasive surgical applications. Given the physical 

* dimensions and material properties of the spring, and the desired 

* deflection where the force will be applied, it will solve for 

* the deflection of the spring on the opposite side of where the force 


* is applied, the force required, and for the angle of deflection. 


* Inputs: Bending stiffness of coil EI, torsional stiffness of the coil 


* JG, number of coils n, radius of the spring R, radius of the spring 
* coil rec, length of spring, and deflection d where the force will be 
* applied. 


* Outputs: deflections d, deflection angle B, and force F. 
Rae ur cE ely wy. G, Roe rey iL. pl. Cc, £1,.aWG,).dE, GOPPF 8B, °Z 
PVTEGER n,m, Kk 


CHARACTER*20 MATL 


PRINT 


PRINT; THIS PROGRAM WILL ALD IN THE DESIGN OF A SPRING FOR USE’ 
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PRINT*,’ IN THE SKELTON FOR MINIMALLY INVASIVE SURGICAL” 


PRINT*,’ APPLICATIONS.” 


pi = 3.1415926 


PRIN ~ 
PRINT*, ‘SELECT AND ENTER THE TYPE OF MATERIAL YOU WISH TO suse. 
PRINT) (1.e. STAINLESS STEEL): ' 
READ* , MATL 
PRINT* 
PRINT*, ‘ENTER MODULUS OF ELASTICITY (E), IN ENGLISH UNITS’ 
PRINT*,"" (i 2e.).28 eo r= 
READ*, E 
PRINT*, ‘ENTER SHEAR MODULUS (G), IN ENGLISH UNITS’ 
PRINE*, *  (isecadte ceo te 
READ*, G 
PR 
PRINT*, ‘ENTER "1" IF YOU HAVE DIMENSIONS IN ENGLISH UNITS (in),' 
PRINT*,*= OR “2™ IF YOU HAVE DIMENSIONS IN SI UNITS (mm) eae 
READ*, k 
PRINT*, ‘ENTER RADIUS OF SPRING: ' 
READ*, R 
PRINT*, ‘ENTER RADIUS OF SPRING COIL: ' 
READ*, rec 
PRINT“, ‘ENTER LENGTH OF SPRING: 
READ*, L 
PRINT*, “ENTER DESIRED DEF LECTICHN a. 
READ*, dF 
IF (k.EQ.2) THEN 
RvR C0 0 Cia) aas Sea) 


me =e BCC 01) tooo 
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ee i O00 Is 837 
Pee Ce Ol) 3037 
ENDIF 
PRINT® 
PRINT*, ‘ENTER NUMBER OF SPRING COILS: ‘ 
READ*, n 
Bein” 


ENE) BNTERe gis 1F PONE FORCE 1s APPLIED,’ 


PRINT*,’ OR "2" IF TWO FORCES ARE APPLIED: ' 
READ*, m 
Poe O2.2o* pi (re**4) 


Tea o pi (re* *4) 


c = n*pi*(R**3) 
Els -E* I 
JG = J*G 


IF (m.EQ.1) THEN 
Ber /(e*(1/(R1) + 3/(dG))) 
dOPPF = c*F*(-1/(EI) + 1/(JG)) 
B = (ATAN((dF - dOPPF) /(2*R)))*(180/pi) 
ELSE 
F = dF/(2*c*(1/(2*EI) + 5/(2*JG))) 
GOPPF = 2*c*F*(-1/(2*EI) + 3/(2*JG)) 
B = (ATAN((dF - dOPPF)/(2*R)))*(180/pi) 


ENDIF 


ERX INT* 

PRINT*, ‘FOR MATERIAL ‘,MATL 

PRINT* 

PRINT*, ‘DEFLECTION AT POINT OF APPLICATION (in): ‘',daF 


emcee, (0200 lL 39.37) 
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PRINT*, ’DEFLECTION AT POINT OF APPLICATION (mm): ‘,dF 
PRINT* 
PRINT*, ‘DEFLECTION AT POINT OPPOSITE APPLICATION (in): ‘,dOPPF 


GOPPF = dOPPE7 (02001739 27) 


PRINT*, ‘DEFLECTION AT POINT OPPOSITE APPLICATION (mm): 3aeEre 
ER oni © 

PRINT*, ‘DEFLECTION ANGLE (degrees): ‘,B 

PRINT*, ’FORCE NECESSARY ((ibsjes) 7 & 

PRINT* 


PRINT*, ‘IF YOU WISH TO BEGIN AGAIN, TYPE ANY NUMBER AND HIT’ 
PRINT*,° RETURN SEU SE ENTER 9399") TO.-OUr la 
READ*, Z 


Tied .NE299) COMO 7> 


END 


The second algorithm developed was for the second model, that of a 
spring bent over directly by an SMA active element. It is a MATLAB 
program where one must input the geometry and material properties of the 
helical spring, as well as the final desired deflection angle of the 
spring. It generates a proportionality matrix whose size is based on the 
number of coils of the spring, and computes the force required to bend the 
spring. It plots force versus angle deflection to obtain a relationship 


between the magnitude of the force required and the amount of deflection. 


&. thesis -m 


%* This program is for the spring design used in the manipulator skeleton 


oe 


for minimally invasive surgical applications. Given the physical 


% dimensions and material properties of the spring, and the desired 
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%$ deflection of the spring, it will solve for the force required to bend 


%$ the spring using an SMA rod inside the spring. 


$ 

teinoucs: - Spring length (L); vadius (RR), coil radius (r), Modulus of 
% Elasticity (E), Shear Modulus (G), number of coils 

$ desired angle of deflection (phi). 


% Outputs: Force required to bend spring (F). 


% 

clear 

% 

IL = input(‘Length of spring in millimeters '); 

R = input(’Radius of spring in millimeters '); 

r = input(‘’Radius of spring coil in millimeters ‘); 


Be=) input (’Modulus of elasticity in GPa ‘); 


G = input(’Shear modulus in GPa '); 


mn = input(‘Number of spring coils (must be even number) '); 
phifinal = input(’Final desired deflection angle in degrees '); 
% 

ee ay 000; 

Ra— h/ 1000; 

pee 2/1000; 

BE = E*leS; 

G = G*le9; 


Paeeinal = phitinal*pi/180; 
% 
Hewes pity 4; 


Jee oe pi ro 4; 


Rom R/(2*E* 1) ; 
Cor 7— &/(2*J*G); 
% 

Pee tos) len 
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Count s=scr 
dradius = 0; 
radius (ijme=- 507 
for Zz =akZ00 


anglephi = L/(2*radius(z)); 


1f anglephi <= phitinal, secon 


Gradius = -0.05*radius(z) ; 


count + 1; end 


radius(z+1) = radius(z) + dradius; 


end 


% 


PC) Baas 0. 
for q = 1-count 
% 
aoe x ayes all 250) 
EOr 9 = 41 ni 
Boat 
end 


end 


EO =-P(q)—k; 
Phitg) = Lyt2se a. 


eV Aelgnl (Vegi igu, 


Clie (2* pit 3) (a*(a 22-4 pie 
(VR — a oly Pe 

C3 0a Fig) 2, 

64 <=) (8"*p172) 7 (aeZ—-¢ "pi 2); 
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e5 


(G4 "pr Sy (ata 2-16" pi 2) }s 


c6 = pi*t(R*2+(P(q)*2/2)); 
gym) (0k 2721) fa = (2*pi*a*R*P(a))/ (a*2=pi-2) 
Lea Zag a) (a) 2, (a val ee 2) 
epee (4 *p1* (a> 2-8 *p1 72) *P(q) *ro)7 (a* (a°2-16"* pi 2) ) 
(4 Ol ach ro), (age 4 «pl 2) 
oem 527 p1°3)/(a*(a>2-16"*pi 72) ); 
el =-sin(a); 
ell = sint(a/2); 
2 = Tom R 2; 
ee = GCl~c3*cl0; 
eae = -e2*c3; 
ela = P(q)*ro*c9*cll; 
ee a=, C/*C10; 
ely j= cl*cl0; 
eo = P(q)*ro*cS; 
perek = 1:n/2 
bore) = 1:k 
Berl, = ls k 
C19 = cos (a*(2*k-1-341)); 
G20 = Sin(a™ (27k-1=-341)); 
2d = cos (a*(i-j)); 
eee = Sin(a*(i=3)); 
eee. = sinta”(k=1+0.5)); 
o24 = sin(a* (k-j+0.5)); 


Mae) ) amon (cls Clo ve eold ten) )s 


Pee ie a emo «Clo cel 0 y= weld e270 = olste23)% 


Poe lear nieMmetels C20 41Cl4 C27 = C15 "e247. 


mie) 


oP? 


A2Z2Z2(i,;j]) = mom* (-cl3*el9 4+ cl4%c7 epi ehee ee 


+ c18*cll*cos((a/2) *(2*k-i-j+41) ) *eos( (a/2) ~*(i-aeee 


Bll(i,j)} = tor* (€6*c2Amaee roel ye 

B12(1,3) = tor*(e6*e22 =sere*c20 2rees- cl ies 

B21(1,3) = tor*(-c6*e22 = Temes ec2Z Ome es ell e724). 

B22(1,J3) = tor*(c6*c21 + cl6*cl19 - c8*cll* (cos (a* (k=1 403s ee 


+ cos(a* (k-3+025))) 4  pi*ro 2); 


Cli(i,)) = el2* (e2*e2ie eile r- 
Cl2(1; 3) = el2* (€276e22-el7 e620). 
G21(1,;)) = vel2* (=e2"e22 ene e200: 


C22(1,3) = €12* (627620 ely eo. 


end 
end 
Git = All+B11+Clil; 
di2 = Al12Z+B12+C12; 
da21 = A21+B21+C21; 
d22 = A22+B22+C22; 


fergiper lei 
Lor] .=22k 
Dit, J) = Dia) eae Gara 
D(i+tk, 7). =o Dia 4k, jee ee il 
D(i, J+K) 2 =i ay eek een calle (er) 
D(it+tk,j+k) = D(i+k,j+k) + d22(i,3); 
end 
end 
end 


B= 32. 
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bor 20> 1n 
BiGre a) oOo Dia 1) 


end 


EQeert— 2 :n/Z 


Snoea( i)? = (n/2 — (1-1)) *a; 

deteania) = (1 -— cos(theta(1))> “dP; 

deltat(i) = (sin(theta(1)) - theta(i))*dP; 
end 


delta = [{deltan’ ;deltat']; 
F = inv(D) *delta; 


Por. = 1:n/2 


fx(i) = F(1i)*cos(theta(i)) - F(n/24+1)*sin(theta(i)); 
fy(i) = F(1i)*sin(theta(i)) + F(n/2+1)*cos(theta(i)); 
end 


Ciel ix’; fy’); 
fe fit Exy; 
Bomee(:,q) = £; 
Prec) = phi(q)*(180/pi) ; 
Sipe=  —-0205"P(q)> 
Edtaj) = P(q) + dP; 

end 

% 

clg 

Pomegn= 1:n/2 
ploerePHl,force();,:)) 
title(’X-Force vs Radius Angle Phi, for Coil #') 
xlabel(’Phi (degrees)') 
ylabel(‘Force (N)’) 


meta thesplot 


=). 


end 
FOr <j =. ny 22a sn 
DpIGE(PHI, force, -)) 
title(’Y-Force vs Radius Angle Phi, for Coil #’) 
xlabel(’Phi (degrees) ’) 
Ylabel(’Force (N)’') 
meta thesplot 


end 


The final element is a set of plots using the previous program. They 
are force versus deflection angle graphs generated by the program using 
the inputs of a 20 mm long spring, with a radius of 5 mm, coil radius of 
0.5 mm, and eight coils. The material chosen was stainless steel with a 
modulus of elasticity of 190 GPa, and shear modulus of 73 GPa. The final 
desired deflection angle was 30°. The eight plots are for the top half of 
the spring. The forces in the bottom half are going to be identical as 
the spring is symmetrical. Of the eight force/deformation plots, four are 


for normal forces, and four are for tangential forces. 
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leoree (IN) 


lorce (IN) 
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